Topological transitions in carbon nanotube networks via nanoscale confinement 
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Efforts aimed at large-scale integration of nanoelectronic devices that exploit the superior elec- 
tronic and mechanical properties of single-walled carbon nanotubes (SWCNTs) remain limited by 
the difficulties associated with manipulation and packaging of individual SWNTs. Alternative ap- 
proaches based on ultra-thin carbon nanotube networks (CNNs) have enjoyed success of late with 
the realization of several scalable device applications. However, precise control over the network 
electronic transport is challenging due to i) an often uncontrollable interplay between network cov- 
erage and its topology and ii) the inherent electrical heterogeneity of the constituent SWNTs. In 
this letter, we use template-assisted fluidic assembly of SWCNT networks to explore the effect 
of geometric confinement on the network topology. Heterogeneous SWCNT networks dip-coated 
onto sub-micron wide ultra-thin polymer channels exhibit a topology that becomes increasingly 
aligned with decreasing channel width and thickness. Experimental scale coarse-grained computa- 
tions of interacting SWCNTs show that the effect is a reflection of an aligned topology that is no 
longer dependent on the network density, which in turn emerges as a robust knob that can induce 
semiconductor-to-metallic transitions in the network response. Our study demonstrates the effec- 
tiveness of directed assembly on channels with varying degrees of confinement as a simple tool to 
tailor the conductance of the otherwise heterogeneous network, opening up the possibility of robust 
large-scale CNN-based devices. 

PACS numbers: 62.25.+g, 85.85+j, 85.35.Kt, 81.07.De, 62.20.Dc 
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Harnessing the highly efficient electronic transport 
in carbon nanotubes continues to be a challenge due 
to difficulties in manipulating and packaging individual 
CNTs PUS]. In instances where devices have been suc- 
cessfully integrated, issues related to device reliability 
and scalability remain largely unresolved. In contrast, 
top-down integration efforts centered around SWCNTs 
networks (CNNs) as electronic material of choice have 
enjoyed considerable success. Several device applications 
have been realized [HE], from gate-modulated thin film 
transistors (TFTs) 7-9 and RF oscillators [TU] to wafer- 
scale integrated circuits on flexible substrates [TTJ [T2] . A 
key limitation of this approach is that current CNT syn- 
thesis routes yield a mixture of metallic and semiconduct- 
ing nanotubes that cannot be easily separated [T3lU5| . 
This intrinsic electrical heterogeneity leads to networks 
with drastically different effective transport characteris- 
tics: as-synthcsizcd networks can be metallic (m-CNNs) 
or semiconducting (s-CNNs). Since the latter response is 
desirable for gate modulation, it is the presence of active 
(metallic) elements that continues to be an issue. 

Purification efforts targeted at increasing the fraction 
of semiconducting nanotubes hold promise yet they in- 
variably compromise the scalability. Alternative strate- 
gies have revolved around suppressing the transport 
through active elements by eliminating the off-state con- 
duction. Efforts to this end involve patterning hetero- 
geneous networks within channels with controlled topol- 



ogy, by i) catalytic growth of CNTs along morpholog- 
ically pattered substrates |16j . or ii) via solution-based 
approaches [I7J [18]. The former rely on chemical va- 
por deposition (CVD) to directly grow CNNs on pre- 
scribed substrates. Network topology can be tuned by 
controlling the density or distribution of the catalytic 
particles and the underlying substrate morphology, re- 
spectively. Random, partially aligned and fully aligned 
networks have been realized via these routes [8]. How- 
ever, the lack of control over the density and distribution 
of catalytic particles and the substrate morphology, espe- 
cially at synthesis temperatures, limits the technique to 
low temperatures and results in network topologies which 
are coupled to the coverage. 

In this article, we use a solution-based approach to 
demonstrate a facile yet scalable route for synthesizing 
CNNs with prescribed transport characteristics. A no- 
table aspect of this approach is the robust control over 
the CNN topology, which we accomplish via template- 
directed fluidic self-assembly of dispersions of functional- 
ized SWCNTs onto suitably patterned channels with pre- 
scribed dimensions (W c x L c ). Plasma treated ~ 100 nm 
silicon dioxide channels are grown thermally on a sil- 
icon wafer and 150 — 500 nm polymethylmethacrylate 
(PMMA) thick photoresists are used to etch nano- to 
micron-scale channels using e-beam lithography. The 
channels are then dip-coated in a dispersion solution 
of commercially acquired SWCNTs (~ 0.25%, mean 
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SWCNT length I = 2.5 ^m, mean SWCNT diameter 
d = 1.25 nm) in dionized water for 5min and then re- 
moved at a rate 100 fiia/s. For these range of parame- 
ters, the receding contact line on the exposed channels 
leads to evaporation-controlled deposition of SWCNTs, 
as illustrated schematically in Fig. flja). For details on 
the synthesis technique, see Ref. [T7J7 

We explore the effect of nanoscale confinement on 
CNNs deposited on patterned channels with varying 
width. Figure [T|b)-(e) shows scanning electron micro- 
graphs (SEM) of the self-assembled SWCNTs for four 
channels widths, W c = 9 /im, 1 /im, 500 nm and 200 nm. 
In each case, the channel length is significantly larger 
(~ 20 /im). As expected, the self-assembly of the func- 
tionalized SWCNTs is restricted to the exposed Si0 2 
channels. Atomic force microscopy (AFM) across the 
channels (not shown) indicates that the assembly is uni- 
form and consists of SWCNT multilayers with thick- 
nesses ranging from a few to tens of SWCNT layers. The 
AFM images also reveal an abrupt (4-6 fold) decrease 
in the thickness as the channels approach nanoscale con- 
finement, defined here as the regime for which SWCNT 
length becomes smaller than the channel width, I < W c . 
Since the SWCNT concentration in the solution and the 
dip-coating withdrawal rate are held constant for all tem- 
plates employed in this study, the thickness change sug- 
gests that the network topology is qualitatively different. 

Direct evidence of the interplay between channel width 
and the CNN topology is obtained from polarized Raman 
spectroscopy on the as-deposited template. The spectra 
are depicted in Fig. [ljf) adjacent to the corresponding 
micrographs. The vertical and horizontal polarization 
intensity curves correspond to the SWNTs aligned along 
the length and width of the channels, respectively. Com- 
parison of the Raman spectra for the different widths 
shows that for confined geometries, I < W c , the peak 
in the vertical intensity curve increases at the expense 
of those in the horizontal and unpolarized peaks - the 
SWCNTs assembly becomes increasingly aligned along 
the channel due to confinement. 

We expect the alignment to also affect electronic trans- 
port along the channels. As a first step, we have ex- 
tracted the device resistivity using 3-probe resistance R c 
(contact + channel) measurements along the channels. 
To calculate the electrical resistivity (= R c W c t/L p ), the 
probe length L p was obtained from SEM images at five 
positions along the channel. Figure [2] shows the variation 
in the resistivity for the four channel widths correspond- 
ing to Fig. [T] As the assembled SWCNT network be- 
comes confined, we see close to two orders of magnitude 
increase in the resistivity consistent with a more aligned 
topology that effectively shields the metallic SWCNTs 
within the heterogeneous network, resulting in transport 
that becomes increasingly controlled by CNT junction 
characteristics (as for example due to CNT bundling). 

The assembled CNNs also vary in thickness and pos- 
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FIG. 1: (a) Schematic illustration of fluidic assembly of SWC- 
NTs onto micropatterned channels consisting of alternating 
pattern lines of Si02 channels (lightly shaded) separated by 
photoresist (dark shaded). The receding contact line formed 
by the SWCNT solution following the dip-coat is exaggerated, 
(b-e) SEM micrographs of the SWCNT-patterned channels 
with four difference widths, (a) W c = 9 /im, (b) 1 /im, (c) 
500 nm and (d) 200 nm. The thickness t of each as-deposited 
channel is the average of fifteen different AFM measurements 
(see AFM images in Supplementary Documents) and indi- 
cated on the corresponding micrograph. The thickness vari- 
ation is due to that in the height of PMMA photoresists. 
Note that the channel width serves as the scale-bar for each 
micrograph. The inset in (a) is a higher magnification mi- 
crograph that shows the detail of the network topology. The 
corresponding polarized Raman spectra are shown below each 
micrograph. Vertical polarization is along the channel. 



3 




100 1000 

Channel width W c (nm) 



10000 



FIG. 2: The effect of channel width on the resistivity of the 
assembled SWCNT networks. While the resistance is the sum 
of the contributions from the probe (contact) and SWCNT 
network, the variation is entirely a network contribution as 
the probe resistance is not width dependent. 



sibly SWCNT density, both of which are expected to 
modify the transport characteristics. To capture the ef- 
fect of these network features in detail, we turn to sys- 
tematic model computations of nanotube assembly on 
experiment-scale channels. SWCNTs coarse-grained as 
rigid rods are employed to gain sufficient statistics on 
the interplay between topology and the nature of electri- 
cal percolation within the network; the effect of channel 
width W c is explored using 2D computations while multi- 
layer, quasi-2D computations are used to study the ef- 
fect of channel thickness. SWCNTs are randomly placed 
onto channels with prescribed dimensions until the de- 
sired network density pnt is obtained. The stochastic 
nature of this procedure allows us to capture the varia- 
tions in network topologies formed due to combination 
of hydrodynamic forces due to the receding contact line 
as well as the evaporation flux that drives the fiuidic as- 
sembly. The average length of the SWCNTs is matched 
to that in the experimental samples, I ~ 2.5 /im and the 
distribution of SWCNT lengths is based on earlier exper- 
imental studies which indicate Weibull statistics [15] . 

The random stick model that has been employed in 
earlier studies [TB] is clearly inadequate in capturing key 
phenomena such as CNT bundling during self-assembly 
(see Fig. [2]). In order to generate realistic topologies, we 
further relax the network by allowing the coarse-grained 
SWCNTs to interact. The inter-CNT interactions, lo- 
calized at the junctions, are obtained by integrating the 
well-known Lennard-Jones (LJ) based description of the 
van Der Waals between graphene surface elements. As 
an example, for fully aligned CNTs, the axially averaged 
inter-CNT interaction energy per length U (R, r) is again 
a 6-12 L J- type potential with constants that are scaled 
by surface integrals which depend entirely on the ratio of 
the inter-tube distance to the CNT radius R/r [2UII2T] . 



Note that this inter-tube interaction is short-ranged and 
negligible for R > y3r, i.e. it is limited to first near- 
est neighbors. For partially aligned CNTs, the van Der 
Waals potential can again be integrated over the surfaces 
of the (pair) of CNTs. The resultant effective inter-CNT 
interactions are angular as they now depend on the de- 
gree of misalignment at the CNT-CNT junction <f>, i.e. 
U = U{R,r,(f>) [Wang and Upmanyu, in preparation]. 
In the case of fully aligned bundles, this inter-tube po- 
tential accurately describes equilibrium inter-tube spac- 
ing, cohesive energy per atom and bulk modulus. The 
interactions serve as inputs for classical dynamical sim- 
ulations aimed at locally relaxing the random network 
with respect to both translational and angular degrees 
of freedom of the individual SWCNTs. Sliding between 
the CNTs at the junctions is ignored. A time step of 1/xs 
is employed and the simulations are performed until the 
interaction energy associated with the network stabilizes. 

Figure [3]ja) shows specific instances of the electrically 
heterogeneous network topologies obtained in three sim- 
ulations with varying widths, W c — 1 /im, 500 nm and 
100 nm. In each case, the SWCNT network is con- 
fined, l/W c < 1 which forces the topology to become 
increasingly aligned with decreasing width. While the 
as-generated random topologies are geometrically aligned 
along the channel, we find that the SWCNT interactions 
always work towards increasing the degree of alignment. 
This is not surprising as the orientation dependence of 
the interaction potential favors a nematic-like phase con- 
sisting of fully aligned nanotubes [22]. For each relaxed 
CNN with a prescribed network density, the electrical 
transport characteristics are extracted by fixing the over- 
all ratio of semiconducting to metallic SWCNTs to the 
theoretical heterogeneous density (3:1). To this end, each 
SWCNT is randomly assigned a metallic or semiconduct- 
ing character and the overall percolation across the CNN 
is measured. The overall percolation can result in i) open 
circuit (OC), ii) semiconducting, or iii) metallic conduc- 
tance across the network. Multiple simulations (~ 100) 
are performed for each channel geometry and density to 
determine the form of percolation. 

Figure [3]^b)-(d) shows the percolation probability 
through monolayered CNNs as a function of the net- 
work density for the three widths shown in Fig. |3ja) . All 
CNNs exhibits two transitions as the network density is 
increased: OC-to-semiconducting at low densities, and 
a semiconducting-to-metallic at high densities. Qualita- 
tively similar transitions are also observed for unrelaxed 
networks that are used as input in random-stick models 
(shown as dotted lines) and comparisons with percola- 
tion in relaxed structures indicates that the enhanced 
alignment driven by the SWCNT interactions shifts the 
transitions to higher network densities. The effect of de- 
creasing width is similar as it also enhances the alignment 
along the channel - for a given network density, the net- 
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FIG. 3: (color online): (a) CNNs observed in three different 2D simulations with varying channel widths, W c = 1 /im, 500 nm 
and 100 nm. In each case, 5 /im length of the 20 fim long channel is shown. For clarity, detailed view of a 1 fim long segment 
of the W = 100 nm channel is also shown. Coloring scheme is based on nature of electrical transport - blue and red indicate 
metallic and semiconducting nanotubes, respectively. The relative number of metallic and semiconducting nanotubes is fixed at 
the theoretical ratio, 1:3. (b-d) Plots of the probability of the nature of electrical percolation through the network as a function 
of network density, Pnt- The solid black, red and blue curves are probability of open circuit, semiconducting and metallic 
conductance across the network, respectively. The dotted lines correspond to simple stick percolation models for randomly 
assembled networks, (e) Contour plot of the probability of semiconducting behavior across the 2D network (shaded red) as a 
function of width (y-axis) and density (x-axis). The length of the channel is fixed at L — 20 /im. (f) Same as in (e) but for 
multi-layer, quasi-2D simulations with varying channel thicknesses, t = 1 — 5 monolayers. The length and width of the channels 
are fixed, L = 5 ^m, W = 200 nm. 



works show a marked reduction in metallic transport as 
the alignment effectively shields the active components. 
This interplay between channel width and network den- 
sity is more clearly seen in Fig. [3^, a contour plot of 
the probability of semiconducting percolation across the 
channel as a function of these two variables. The network 
density ranges for semiconduction increases substantially 
as the width is decreased below W c ~ 300 nm. Note that 
the probability of an open circuit also increases with de- 
creasing width, but the extent of this effect is smaller 
than the enhancement in semiconduction. 

The thickness of the CNNs permits control over out- 
of-plane confinement of the SWCNTs and we explore 
this effect by performing percolation studies on relaxed 
multilayered networks with fixed width and length. For 
computational efficiency, we have chosen networks with 
smaller lengths (L c — 5fim) and channel width (W c = 
200 nm). Note that the network is already confined at 
this width. Figure [3] shows the contour plot associated 
with probability of semiconduction along the channel as 
a function of number of monolayers (ML), t = 1 — 5 ML. 
Our results show that thin CNNs have a higher probabil- 
ity for semiconduction and the probability decreases (at 
the expense of metallic behavior) rapidly within the first 
few monolayers. While we have not studied thicknesses 
greater than t = 5ML, the contour plot shows that the 
marginal decrease should be much smaller. 

In summary, this article demonstrates a scalable 
solution-based route for large-scale integrated circuits 
based on directed assembly of CNNs. Our results show 
the effectiveness of geometrical confinement induced by 



controlling the channel geometries, mainly their width 
and depth. Small, nanoscale geometries suppress the ef- 
fect of active elements inherent in heterogeneous SWCNT 
solutions, and the challenge of large-scale integration 
of CNN-based devices reduces to scalable, high-fidelity 
lithography of assembly patterns. 
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